arXiv:1507.05750vl [quant-ph] 21 Jul 2015 


A chiral route to spontaneous entanglement generation 

Carlos Gonzalez-Ballestero/’ * Alejandro Gonzalez-Tudela,^ Francisco J. Garcia-Vidal/’^ and Esteban Moreno^ 

^Departamento de Ftsica Teorica de la Materia Condensada and Condensed Matter Physics Center (IFIMAC), 
Universidad Autonoma de Madrid, E-28049 Madrid, Spain 
^Max-Planck-Institut fiir Quantenoptik Hans-Kopfermann-Str. 1. 85748 Garching, Germany 
^ Donostia International Physics Center (DIPC), E-20018 Donostia/San Sebastian, Spain 

We study the generation of spontaneous entanglement between two qubits chirally coupled to a 
waveguide. The maximum achievable concurrence is demonstrated to increase by a factor of 4/e ~ 1.5 
as compared to the non-chiral coupling situation. The proposed entanglement scheme is shown to be 
robust against variation of the qubit properties such as detuning and separation, which are critical 
in the non-chiral case. This result relaxes the restrictive requirements of the non-chiral situation, 
paving the way towards a realistic implementation. Our results demonstrate the potential of chiral 
waveguides for quantum entanglement protocols. 
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Efficient quantum circuits are a very important ingre¬ 
dient for the development of quantum computing^. Usual 
implementations of these devices require platforms where 
information, usually in the form of photons, can be easily 
introduced and extracted^. Several systems based on 
waveguides have been proposed for quantum circuitry, 
from superconducting stripes^ to dielectric^’^, photonic 
crystal®’^ or plasmonic waveguides®. In this context, the 
interaction between the guided photons and quantum 
emitters is critical for various processes such as the cre¬ 
ation of entangled states between the qubits. Spontaneous 
entanglement generation in waveguide setups has already 
been predicted®’Many other interesting phenomena 
like mesoscopic entanglement^^, long-distance quantum 
beats^®, or the formation of sub- and superradiant states^®, 
show that waveguides are excellent platforms for quantum 
information processing. 

Recently, systems of emitters chirally coupled to waveg¬ 
uides have attracted a lot of attention both theoretically 
and experimentally^^^^®. In these configurations, an ade¬ 
quate engineering of the waveguide can be used to break 
the emission symmetry of the qubits, channelling the 
emitted photons preferentially into one of the two direc¬ 
tions of the waveguide. Chirality in emitter-waveguide 
coupling is a general effect associated to the so-called spin- 
orbit interaction of light^^. Besides theoretical studies, 
many experiments have reported chiral emission with a 
very large degree of directionality, from nanoparticles and 
atomic ensembles in dielectric waveguides^^’^® to quantum 
dots in nanobeams^® and photonic crystals®®^®®. Among 
these setups, the latter turn out to be especially promising 
systems as they combine large directionalities of around 
90% with high emitter-waveguide coupling fractions (up 
to 98%). Consequently, they have been proposed as ideal 
platforms for implementation of quantum logical gates®®. 

In this Rapid Communication, we tackle the problem of 
spontaneous entanglement generation between two qubits 
chirally coupled to a waveguide. In the first part of this 
work, we present an analytical solution to the Master 
Equation describing the evolution of the system state, 
showing how chirality allows for an enhancement up to 



FIG. 1. System under study. Two qubits of frequency cjq 
and separated by a distance d are placed in the vicinities of 
a waveguide. The energies jja quantify the chiral coupling 
of qubit j to the photonic propagating mode a (= L, R). In 
the same fashion, the decay rate into the 3D environment and 
other lossy modes is described by a coupling constant Fj. 

~ 50% in the maximum generated entanglement as com¬ 
pared to the non-chiral case. In the second part of the 
paper, we present a more complete formalism in which 
non-Markovian effects are explicitly accounted for by fully 
diagonalizing the Hamiltonian in the single-excitation 
subspace. We use this formalism to demonstrate the ro¬ 
bustness of the entanglement generation scheme against 
the detuning between the qubits, the total coupling rate, 
and the qubit-qubit separation. 

The system under study is depicted in Fig. 1. Two emit¬ 
ters I and 2, modelled as two-level systems of frequency 
Wo, are coupled to the propagating photonic modes of 
a waveguide. The emitters are separated by a distance 
d = X 2 — Xi, and coupled to the right and left propagat¬ 
ing photons through the energy constants jja and 7jL 
(j = 1,2), respectively. In a chiral coupling scheme such 
as the one analyzed in this work, these constants are dif¬ 
ferent {jjn 7 ^ Tjl)- Other deexcitation processes into free 
space or additional lossy modes are taken into account 
through the decay rates T^. The three coupling constants 
of each qubit are used to define a usual figure of merit in 
waveguide systems, namely the coupling fraction or beta 
factor, given by /3j = (Tja -b Jjl)/ {IjR + + F^). 

Our aim is to analyze spontaneous entanglement gen¬ 
eration when qubit 1 (the qubit on the left) is initially 
excited. We start by solving the system dynamics, ob- 
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taining the time evolution of the reduced density matrix 
of the two-qubit subsystem, p. In the first part of this 
work, we will follow a usual approach undertaken in quan¬ 
tum optics, in which the problem is simplified by tracing 
out the photonic degrees of freedom under the so-called 
Markov approximation^^. In this situation, the only dy¬ 
namical variable is the density matrix p, whose evolution 
is governed by the following general Master Equation; 

p=-i [H, p]+ [p] + 

i=i.2 (I) 

[p]+V72Z,7ii'^o-i,(T2 [p]) 

where we define jj = (jjR + 'jjL) /2. The generalized 
Lindblad superoperators, £ 0 -^ , are employed to describe 

incoherent processes, in this case a waveguide-mediated 
interaction. They are expressed as 

- e^^-^-nala,p]), (2) 

where Dab = ja^a — a^bl/Ao, Aq = 2'KVg/ujQ is the emission 
wavelength of the qubits, and Vg is the group velocity 
of the guided photons. The bare Hamiltonian of the 
system is given by = wo(cr|cri a\a 2 ), where Oj is the 
annihilation operator of qubit j. Note that Eq. (1) is 
particular of chiral configurations^®, and is reduced to its 
more common form in the case = "/jL- In order to 
focus on the fundamental aspects of the chiral system, we 
will first particularize our study to the lossless case (i.e., 
f3j = 1), including the losses in the second part of the 
work. 

Combining Eq. (1) with the particular initial conditions 
of our problem, p(0) = CT||0)(0|cri, and expressing p in the 
usual basis {|0),|I) = cr||0),|2) = crj|0), |3) = o-|crj|0)}, 
the only non-zero elements of the density matrix are the 
populations poo, Pii, P22, and the coherence pi2. Three of 
these quantities are coupled through the following system 
of differential equations, 

Pii = 

— 271P11 — v^ 7 il 72 l , 

P22 = 

— 272 P 22 — -f , 

P12 = —(71 + 72)^12 

~V7TRT^Plie *27r!i_ 

with the normalized distance d = d/Xo. Once these 
equations are solved, we can compute the qubit-qubit 
entanglement which we quantify with the Wootters con¬ 
currence (7^^, a widely used measure thanks to its simple 
calculation and its intuitive bounds. Indeed, the concur¬ 
rence ranges from 0 for non entangled states, to 1 for 
maximally entangled configurations. A straightforward 
calculation demonstrates that in this case, C = 2|pi2|. 



d/ Aq 


FIG. 2. Time evolution of the concurrence, Eq. ( 6 ), for differ¬ 
ent qubit-qubit separations d, in the non-chiral (a) and chiral 
(b) cases, c) Dependence of Cmax on the separation between 
the qubits, d. The condition Ai = A 2 is chosen because it op¬ 
timizes the maximum achievable concurrence Cmax (see main 
text). In all three panels the ideal case Pj =1 is considered. 


An analytical solution of Eq. (1) can be obtained when 
the two qubits are equally coupled, i.e., 71 = 72 = 7 - The 
following expression for the concurrence is obtained: 


C^(t) - i I + + 47t 


( 6 ) 


X sm 


2q"ft sin(27rd) 


-I- sinh 


2q'yt cos(27rd) 


where g = (l — (l — and we have intro¬ 

duced the directionality of qubit j as the adimensional 
ratio Aj = {-fjR - -fjL) / iljR + Ijl)- 

In order to have a clear view of the advantages brought 
about by chirality, let us first analyze the non-chiral case, 
Aj = 0. The time evolution of the concurrence in this 
situation is displayed in Eig. 2a. Regardless of the sepa¬ 
ration between the qubits, the concurrence is 0 for t = 0 
since the initial state is separable, and grows up to a 
maximum value as, due to the exchange of population 
through the waveguide modes, the state of the qubits 
becomes entangled. Eor t —>■ 00 , the behavior of the 
concurrence depends on the separation d. Eor almost 
all separations, the population abandons the vicinities 
of the emitters in the form of propagating photons and, 
as a consequence, the concurrence decays to zero. How¬ 
ever, when the qubits are identical and for very specific 
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separations, 2d = 0, 1 , 2,the behavior is different, as 
the concurrence not only reaches its maximum achievable 
value, C = 0.5, but also has an infinite lifetime. This 
peculiar time evolution is caused by the appearance of a 
Fabry-Perot-like resonance between the qubits, where a 
photon of energy wg can be trapped forming a standing 
wave^®. The presence of this localized photon is linked 
to the fact that the transmittance of one qubit is strictly 
zero for an incoming resonant photon^®. Note that, in a 
realistic case (/3 < 1) the concurrence always decays with 
time^°. 

When the coupling is chiral, on the other hand, a 
straightforward calculation demonstrates that the above- 
mentioned standing wave does not appear^^, as the chi¬ 
rality effectively couples the right- and left- propagating 
modes and, consequently, the single qubit transmittance 
never vanishes. As figure Fig. 2b shows, this results in a 
different time evolution in which two important features 
arise, namely the weak dependence on the separation d, 
and the achievement of much larger values for the con¬ 
currence as compared to the non-chiral situation. The 
maximum concurrence achieved during the time evolution, 
Cmax, is displayed in Fig. 2c as a function of d. Clearly, 
the maximally chiral configuration Aj = 1 optimizes 
the entanglement generation scheme, as Cmax reaches a 
maximum value which, additionally, is independent on d. 
Chirality thus represents an advantage towards a realistic 
implementation, as it can overcome the critical depen¬ 
dence on the separation in non-chiral configurations. 

In spite of the weak dependence with the qubit-qubit 
separation, there exist nevertheless specific separations 
which are optimum for the entanglement generation. As 
Fig. 2c shows, the concurrence Cmax is maximized when 
the distance between the qubits is 2d = 0,1, 2... For these 
separations, an analytical expression can be extracted 
from Eq. (6), 

f^2 _ (1 + Ai)(l -I- A 2 ) / 1 - 

l-g2 \l + q) ' ^ ’ 

To understand the dependence of Cmax with the direction¬ 
alities, which is shown in Fig. 3, it is important to bear in 
mind that initially only qubit 1 is excited. For A^- = — 1 
the concurrence is strictly zero during all the time evo¬ 
lution, as the qubit j is coupled only to left-propagating 
modes and thus is not able to interact with its partner. 
On the other hand, any values of the directionalities in 
the region Ai, A2 > 0 result in an enhancement of Cmax 
with respect to the non-chiral case. Moreover, when both 
qubits are maximally coupled to right-propagating modes, 
i.e., Ai = A 2 = 1, the maximum entanglement rises up 
to a very large value, which can be extracted from Eq. (7) 
as lim,j_,,o Cmax = 2/e ~ 0.73. This is a significant result, 
as it shows that the maximum achievable concurrence can 
be enhanced by a factor of ~ 50% with respect to the 
non-chiral coupling scheme. 

We have seen that playing upon chirality it is possible 
both to reach higher values of the maximum concurrence. 


max 
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FIG. 3. Maximum achievable concurrence in the ideal case 
{j3 = 1 ), versus the directionalities of each qubit, Ai, A 2 , for 
the left qubit initially excited, and a separation d = Aq. The 
blue line displays the non-chiral value Cmax = 0.5. 

Cmax, and to reduce the sensitivity to the qubits separa¬ 
tion, d. However, for large separations or coupling rates, 
'yd^ Vg, it has been reported that non-Markovian effects 
arise in waveguide systems^^. Such effects introduce ad¬ 
ditional retardation that cannot be described with the 
currently used master equation formalism. Hence, we now 
reformulate the problem to properly assess the robustness 
of the proposed protocol. We employ a more complete 
approach by diagonalizing the full Hamiltonian of the 
system in the single excitation subspace, as detailed in 
the supplemental materiaP^. Unless stated otherwise, we 
will assume equally coupled qubits, i.e., = 72 ^ and 

7 iL = l 2 L and large directionalities (Ai = A 2 = 0.90) in 
order to stay close to the optimum configuration. We also 
include explicitly the lossy modes in this part of our work, 
by fixing the decay rates Fj such that jdj = 0.98. Note 
that the chosen values for both the directionalities and 
the beta factors have been experimentally reported^^. 

We employ the above mentioned formalism to explore 
the robustness of the scheme against variation of three 
parameters, namely the detuning between the frequen¬ 
cies of the emitters, the total coupling to the waveguide 
modes, 7 , and the separation d. In Fig. 4a, the effect of 
detuning between the qubits is studied. For the corre¬ 
sponding calculations, the frequencies of both emitters are 
modified according to wi = wg + <5/2, and W 2 = wg — (5/2. 
Physically, we are shifting away the emission spectrum 
of qubit I from the absorption spectrum of qubit 2, keep¬ 
ing their linewidth constant. It is known that a strong 
overlap between both the spectral distribution of the pho¬ 
ton emitted by qubit I, and the absorption spectrum of 
qubit 2, is key for entanglement generation^®. Hence, 
the concurrence naturally decreases for large values of S. 
The detuning relative to cug seems to be the most critical 
parameter as a change of ~ 0.5% in the frequency of the 
qubits is enough to reduce Umax below 0.5. However, the 
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robustness against detuning is considerably large with 
respect to the qubits linewidth, 7 . Indeed, whereas in 
the non-chiral case the concurrence is independent on the 
detuning for 5 ^ O. 27 , for chiral couplings this range is in¬ 
creased by a factor of ~ 5. Additionally, in the chiral case, 
concurrences of Cmax = 0.5 are possible for detunings as 
large as <5 ~ 57 . This is a crucial advantage with respect 
to non-chiral systems, especially for quantum emitters 
with a very narrow linewidth such as quantum dots. 

The variation of Cmax with the total qubit-waveguide 
coupling is displayed in Fig. 4b. For low values of 7 , 
the concurrence is close to its theoretical maximum 
(Cmax ~ 0.7) due to the large chosen directionalities. 
When 7 is increased, the photonic wavepacket emitted by 
qubit 1 becomes narrower in space^* and, eventually, its 
width becomes smaller or comparable to the qubit-qubit 
separation d. As a result, the qubit 1 significantly decays 
before the photon reaches qubit 2 and, at any given time, 
at least one of the emitters is fairly depopulated. Thus, 
the concurrence C = 2 |pi 2 | = 2y/piip22 decreases for 
large values of 7 . This effect has been studied in detail 
in non-chiral conhgurations^^. Note that, nevertheless, 
even for couplings as large as 7 « O.Iwq, the concurrence 
remains above 0.6. This result shows that chirality allows 
for a high level of concurrence not only in the optical 
regime, but also in systems where much larger couplings 
arise such as superconducting stripes. 

Finally, the variation of the concurrence with the qubit- 
qubit separation d is shown in Fig. 4c. While in the 
non-chiral case a maximum concurrence of 0.5 was ob¬ 
tained only for particular values of d, for chiral couplings 
the entanglement generation scheme is shown to be ro¬ 
bust for a wide range of separations. For large distances 
there is a decay in the concurrence, which responds to 
the same mechanism discussed above. In this case, the 
spatial extension of the emitted photon is made smaller 
than the separation d by directly increasing the qubit sep¬ 
aration instead of the coupling 7 . Interestingly, for qubit- 
waveguide couplings in the optical regime (7 < , 

the entanglement generation scheme is extremely robust 
with respect to the distance d, allowing for concurrences 
above 0.6 for very large separations, e.g. around 60 ^m 
at Wo 2 eV. The separation between the qubits is thus 
not a critical parameter anymore, allowing for a much 
easier implementation of this entanglement protocol. 

In conclusion, the phenomenon of spontaneous genera¬ 
tion of entanglement between two qubits chirally coupled 
to a waveguide has been analyzed in detail. We show 
that even the slightest directionalities in the couplings 
may improve the maximum achievable entanglement, as 
compared to non-chiral systems. Moreover, we identify 
the optimal directionalities and demonstrate a very signif¬ 
icant enhancement of the maximum entanglement. This 
entangling scheme displays a fairly weak dependence on 
the relevant parameters, which highlights the robustness 
of the protocol. In particular, when compared to the non- 
chiral case, the influence of the qubit-qubit separation is 
reduced, which constitutes an important advantage for the 
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FIG. 4. Maximum concurrence Cmax as a function of the 
various system parameters (A^ = 0.9, Pj = 0.98). a) Effect of 
the detuning between the transition frequencies of the qubits, 
S. b) Dependence on the total coupling to the waveguide, 7 . 
c) Dependence on the qubit-qubit separation d, for different 
values of 7 . All the coupling constants are normalized to wo. 


feasibility of an experimental implementation. In addition 
to the interest of our specific results, our work positions 
entanglement as one further reason for the exploration of 
chiral waveguide QED. 
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A chiral route to spontaneous entanglement generation: Supplemental Material 


DIAGONALIZATION OF THE HAMILTONIAN. 


We make use of the real-space formalism which is 
usually employed in low-excitation problems in waveg¬ 
uide QED^. The Hamiltonian of the system is H — 
Hqb + Huig -I- Hj , where the two first terms correspond 
to the energy of the qubits and the waveguide, respec¬ 
tively, and Hj is the interaction term. They are given by 
{h=l): 


Hqb = 


Hwg — ^^5 I dx 


cl{x)dj:CLix) - c\{x)dxCii{x) 


( 1 ) 

( 2 ) 


Hi = YY jdx5{x-Xj)[Vjaci{x)ag+h.c]. (3) 

A — 1 ^_D T ^ 


j—la.—R,L' 


In the above equations, 0^- is the transition frequency of 
qubit j, and Vg is the group velocity of the guided modes, 
whose dispersion is considered linear. The constants Vja, 
assumed real for simplicity, are related to the coupling 
rates in the main text through jja = V^^/vg. The oper¬ 
ators tJj and Ca(x) destroy an excitation in qubit j and 
a a— propagating photon at position x, respectively. Fi¬ 
nally, Xj = zizdl2 is the position of the emitter j along 
the waveguide. Note that the losses Fj are not accounted 
for in Eq. (1), as we introduce them a posteriori. 

In order to diagonalize the system Hamiltonian in the 
single excitation subspace, we assume an eigenstate of the 
form 


k) = ( dx(l>a{x)cl{x) j |0), (4) 


and solve the time-independent Schrodinger equation, 
H|e) = e|e), to obtain the coefficients {aj,(j)a{x)}. Hol¬ 
lowing the usual approach, we make a plane wave Ansatz 
for the wavefunctions (paix), i.e. 


f A for a: < —d/2 
B for -d/2 <x<d/2, 
i C for d/2 < X 


(5) 


( D for a: < —d/2 
E for -d/2 <x< d/2 
\ F for d/2 < X 


( 6 ) 


which reduces the problem to an algebraic system of 
equations^. For each energy e, we can naturally find two 
linearly independent eigenstates, corresponding to the 
scattering of photons coming from either a: = — oo or 
X = oo. These states, labeled |e+) and |e_), are obtained 
by setting {A = 1, E = 0} and {F = 1, A = 0} in Eqs 5-6, 
respectively. The eigenstates |e±), whose degeneracy is 
inherited from that of the right/left propagating modes of 
a bare waveguide, have been called scattering eigenstates 
in the literature. 

The scattering solutions of Schrodinger equation have 
been reported to not always form a complete basis, as 
localized resonances may arise^. We thus need to check 
our equations for orthogonal solutions, i.e. A = E = 0. 
A simple calculation shows that the necessary conditions 
for a localized eigenstate are 71 ^ = 71 L and 72 ^: = 72 L- 
We conclude that, in the chiral case, the basis {|e±)} is 
complete, as no localized resonances appear. 

With the eigenstates at hand, we are able to introduce 
the qubit losses, Fj. As localized eigenstates are not 
present in the chiral case, the effect of the lossy modes 
is easily accounted for. Indeed, it has been shown in 
ref.^ that the effect of the external decay of the qubits 
over the scattering eigenstates can be fully reproduced by 
adding an imaginary part to the frequency of the emitters, 
ilj ilj - iFj/2. 

The final step for obtaining the system dynamics is 
to construct the time evolution operator, U{t). If the 
coupling is chiral, this task is not as straightforward as 
in previous works, as the two scattering branches are 
not orthogonal. The operator U{t) takes now a more 
complicated form: 

*e-“l«)(S--),,fel. (7) 

s .,,=±7-“ 

where we define the overlap matrix S as 


S,g = lim ( 8 ) 

L —^00 L 

Using Eq. (7) we can numerically obtain the evolution 
of any initial state through |'0(t)) = 17(t)|V'(0)). Finally, 
it can be checked that C/( 0 ) = 1 , which certifies the 
completeness of the basis {|e±)}. 
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